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Abstract. We present a catalogue and atlas of low-resolution spectra of a well defined sample of 341 objects in 
the FORS Deep Field. All spectra were obtained with the FORS instruments at the ESO VLT with essentially 
the same spectroscopic set-up. The observed extragalactic objects cover the redshift range 0.1 to 5.0. 98 objects 
are starburst galaxies and QSOs at z > 2. Using this data set we investigated the evolution of the characteristic 
spectral properties of bright starburst galaxies and their mutual relations as a function of redshift. Significant 
evolutionary effects were found for redshifts 2 < z < 4. Most conspicuous are the increase of the average CIV 
absorption strength, of the dust reddening, and of the intrinsic UV luminosity, and the decrease of the average 
Lya emission strength with decreasing redshift. In part the observed evolutionary effects can be attributed to an 
increase of the metallicity of the galaxies with cosmic age. Moreover, the increase of the total star-formation rates 
and the stronger obscuration of the starburst cores by dusty gas clouds suggest the occurrence of more massive 
starbursts at later cosmic epochs. 
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1. Introduction 

The advent of the lOm-class telescopes has allowed di- 
rect access to the early stages of galaxy evolution using 
spectroscopic methods. Observing with Keck, Steidel et al. 
p996allT9^6Bjl and Lowenthal et al. l(TM7|> succeeded first 
in acquiring spectra of high-redshift galaxies at z ~ 3. Part 
of the first objects were selected from the Hubble Deep 
Field North (Williams et al. 119961 review by Ferguson 
et al. I2000[) . All candidates were identified using a two- 
colour selection method based on the Lyman-limit break. 
In the following years many more 'Lyman-break galax- 
ies' were observed at redshifts z ~ 3 (e.g. Cristiani et 
al. 120001 Vanzella et al. 120021 Steidel et al. 12003) . z ~ 4 
(Steidel et al. HOOOjl and higher (Lehnert & Bremer 12003) . 
An extended review of this topic has been presented by 
Giavalisco l|2002|l . 

In addition to the two-colour diagram method, which 
is biased towards objects with high star-formation rates 
and/or strong intergalactic absorption, deep multi-band 
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photometric surveys resulting in 'photometric redshifts' 
(Connolly et al. ITW71 Fernandez-Soto et al. ITrM) have 
also been used to identify high-redshift galaxies. In this 
way, e.g., candidates with 1.5 < z < 2.5 (Savaglio et 
al. 120041 Daddi et al. 12004) and z > 5 (Spinrad et al. 
1998; Weymann et al. I1998J) could be successfully iden- 
tified. Finally, objects with strong Lya emission and a 
weak continuum can be identified using narrow-band fil- 
ters (Hu & McMahon HOOBI Cowie & Hu HOOSI Hu et al. 
ITOOOI Rhoads et al. 120031 Maier et al. 12003) . 

Because of the faintness of high-redshift galaxies, de- 
tailed spectral studies have so far been restricted mainly 
to composite spectra (Lowenthal et al. 119971 Steidel et al. 
2001; Shapley et al. 2003, S03 hereafter) or to galaxies am- 
plified by gravitational lensing (Pettini et al. 120001 [5002 
Mehlert et al. l200Tl Frye et al.|2002 Hu et al. l2002|) . 

The various spectroscopic studies of Lyman-break 
galaxies have shown that these objects are vigorously star- 
forming galaxies having large-scale outflows of neutral gas. 
Typical outflow velocities of a few lOOkm/s were derived 
from the blue-shift of the interstellar absorption lines (e.g. 
Pettini et al. 2000] Adelberger et al. 12003) . Investigations 
of the shape of the Lya line and its relation to other spec- 
tral properties fsee lS03|l revealed a complex dependence of 
the Lya emission strength on the (neutral) gas and dust 
distribution and kinematics in these galaxies. Pettini et al. 
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l|2000l I2002[l using weak low-ionisation lines and Mehlert 
et al. H2002jl using the high-ionisation stellar-wind blend 
C IV AA 1548, 1550 find evidence for a chemical evolution 
with cosmic age of the bright starburst galaxies. Further 
indications for evolutionary effects come from photomet- 
ric data on the luminosity, size and mass of high-redshift 
galaxies (e.g. Lowenthal et al. 119971 Madau et al. 119981 
Steidel et al. USSS! Shapley et al. ESUUTl Ferguson et al. 
2004 Idzi et al. 2004), suggesting the formation of more 
massive galaxies at later cosmic epochs. 

In the present paper we describe a deep spectroscopic 
survey using high-quality low-resolution spectra of galax- 
ies covering the redshift range < z < 5. This survey al- 
lows a detailed analysis of the dependence of basic galaxy 
properties on cosmic age. We primarily focus on the anal- 
ysis of evolutionary effects in the redshift range 2 < z < 4, 
where we expected the most significant results. A detailed 
discussion of other redshift ranges will be presented else- 
where. 

This survey is part of the FORS Deep Field (FDF) 
project (Appenzeller et al. 2000; Hcidt et al. 2003b) which 
has been carried out at the ESO Very Large Telescope 
(VLT) as a guaranteed time programme using the two 
FORS instruments. The FDF is located close to the South 
Galactic Pole and has a size of about 7' x 7'. The photom- 
etry is based on deep images in nine filter bands from U to 
K (Heidt et al. 12003 bl Gabasch et al. l2tMjl . which allows 
an efficient selection of candidates for spectroscopy using 
photometric redshifts. The 50% completeness limits in the 
Vega system are 27.7, 26.9, 26.7, 26.4 mag in B, g, R, and 
/, respectively. 

In the following we describe the FDF spectroscopic 
survey, discussing the sample selection (Sect. , the ob- 
servations (Sect.0), the data reduction (Sect.QJ and the 
derivation of redshifts and object types (Sect.EJ). The cat- 
alogue of the FDF spectroscopic sample (available in elec- 
tronic form only) is described in Sect. El Basic characteris- 
tics of the spectra and the spectroscopic redshift distribu- 
tion are outlined in Sect. and |HJ respectively. In Sect. [5] 
we analyse the FDF high-redshift sample with regard to 
evolutionary effects. Implications from the observational 
results are discussed in Sect. 1101 

Throughout the paper H = 67kms~ 1 Mpc -1 , Q\ = 
0.7 and f^M = 0.3 are adopted. 

2. The spectroscopic sample 

For the present investigation we selected a sample of galax- 
ies in the FDF covering a large range of photometric 
redshifts derived from the FDF photometric programme 
(Heidt et al. ;2003b). The photometric redshifts were cal- 
culated by fitting semi-empirical spectral energy distri- 
butions (SEDs) to the filter-band fluxes (Bender et al. 
I3UUT1 I2UU1 Gabasch et al. E?tM|l . The templates were 
constructed using stellar population models of Maraston 
( 199SJ) with variable reddening according to Calzetti et 
al. I)1994|l and broad-band SEDs produced from photo- 
metric data of galaxies in the Hubble Deep Field North 
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Fig. 1. / magnitudes as a function of the photometric red- 
shifts of the primary (circles) and secondary extragalac- 
tic targets (crosses) of the FDF spectroscopic survey. The 
bright object with z p hot = 3.1 is the quasar Q 0130-260 
(Zspec = 3.365). 

(HDF-N, Williams et al. I1996| with known spectroscopic 
redshifts. From a comparison of photometric and spectro- 
scopic redshifts carried out in the HDF-N and the FDF 
we estimate for the present FDF spectroscopic sample a 
fraction of misidentifications of ~ 3 % only. The rms error 
of Zphot — z S poc for the FDF objects with spectroscopically- 
confirmed photometric redshifts amounts to 0.13. 

To obtain sufficiently large subsamples of candidates 
for the different redshift intervals, apparent brightness 
limits for the object selection were chosen depending on 
the redshift. Fig.^shows the / magnitudes (Vega system) 
of the spectroscopically-observed objects as a function of 
photometric redshift. For the redshift range 2.0 to 4.0 the 
selection limit was normally 24.5 mag in I. For / < 24 mag 
the fraction of the spectroscopically-observed photometric 
candidates was about 50%. For lower redshifts the same 
level of completeness was achieved for brighter objects 
(/ < 22.5 mag for 1.0 < z pho t < 1-5). For z pho t > 4.0 
one third of the objects sufficiently bright for spectroscopy 
(J < 25 mag) was included. 

The majority of the targets selected are intrinsically 
bright objects (L > L*) at redshifts between 1 and 5. In 
addition to these primary targets the spectroscopic sam- 
ple contains objects which were observed serendipitously 
as their images coincided by chance with the slit positions 
of the primary targets. Two or three additional objects on 
a slitlet were not unusual. These secondary objects con- 
stitute a random sample constrained only by the appar- 
ent brightness. Hence, faint intermediate-redshift galaxies 
dominate this sample, as indicated in Fig. ^ 

A few high-redshift galaxy candidates were selected 
which deviated from the selection criteria for z > 2 pri- 
mary targets. Firstly, three Lya emitters, which were ini- 
tially observed as secondary targets, were converted to 
primary targets. Secondly, we selected three faint high- 
redshift candidates located in the vicinity of z p hot > 2 
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primary objects with / < 24.5 mag. Thirdly, we included 
three faint objects located in the vicinity of the bright 
FDF quasar Q 0130-260 (see Fig. EjJ. These three can- 
didates had been selected as likely Lya emitters from 
narrow-band images sensitive to Lya at the quasar red- 
shift (see Heidt et al. I2003a|l . Furthermore, three objects 
were selected on the basis of their photometric redshifts 
2 p hot ~ 3 and a colour excess g — R> 0.5 mag, indicating 
strong Lya emission. In principle, these additionally se- 
lected objects increase the fraction of Lya bright galaxies 
in the sample. However, the sample based on the standard 
brightness criterion is not significantly affected, since the 
continuum flux of the Lya emission galaxies is generally 
low. Only two of the additional galaxies are brighter than 
I = 24.5 mag. Consequently, the high-redshift sample can 
be considered as a representative sample for the absolute 
brightest objects with z > 2. 

3. Observations 

All spectroscopic observations of the FORS Deep Field 
were carried out using the two FORS instruments at the 
ESO VLT at Paranal, Chile, mainly between September 
2000 and September 2002. The total integration time was 
63. 7h. 

All spectra were taken using the multi-object spec- 
troscopy modes MOS and MXU of FORS (see FORS 
Manual). The MOS masks consist of 19 individually mov- 
able and adjustable slitlets. The MXU masks are manufac- 
tured using a laser-cutting device. Since the MXU was not 
available for the first campaigns and only at FORS 2, al- 
most 90 % of the exposures were obtained using the MOS 
mode. 

To assure a homogeneous sample the observations were 
carried out using a standardised set-up. The slit width 
was uniformly set to 1", corresponding approximately to 
the size of typical high-redshift galaxies under average at- 
mospheric conditions at Paranal. For all observations dis- 
cussed here, the low-resolution grism 1501 was used. This 
grism covers the full range of the CCD sensitivity (3300 - 
10000 A) with a relatively high efficiency, reaching its max- 
imum at around 5000 A. The grism was employed without 
order separation filters in order to exploit the whole wave- 
length range and to make most efficient use of the observ- 
ing time. Since most galaxies were red with no or little blue 
flux, second order contamination in the red was normally 
negligible. In the case of blue objects with significant sec- 
ond order contamination only the uncontaminated part 
of the spectrum was used. The set-up resulted in a mea- 
sured spectral resolution element of w 23 A (FWHM) for 
the calibration lines and sa 18 A for point sources under 
good seeing conditions. 

Our high-redshift galaxies typically had an apparent 
brightness of I > 23 mag. On the other hand, signal-to- 
noise ratios of 10 and more are needed for quantitative 
spectroscopic studies. Therefore, integration times up to 
10 h were scheduled to achieve the desired S/N. The obser- 
vations of the individual objects were usually distributed 
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Fig. 2. Continuum signal-to-noise ratio per resolution cl- 
ement of the FDF spectra (in the filter band B, g, R or I 
with the best average S /N) as a function of the I magni- 
tude. Circles refer to primary targets, crosses to secondary 
targets. 

over several mask set-ups to optimise the observing time 
allocation and to reduce systematic errors caused by CCD 
defects and object locations near slit edges. The individual 
exposures had integration times between 30 and 48 min. 

4. Data reduction 

The data were reduced using MIDAS 1 routines optimised 
for the reduction of FORS data, mainly based on the 
MIDAS context MOS. 

For the bias subtraction for each observing run a mas- 
ter bias was generated by calculating the normalised me- 
dian of a minimum of five single bias frames and sub- 
tracted from the raw frames after scaling the bias level 
with the overscan. 

The pixel-to-pixel sensitivity variations were corrected 
by means of internal screen flat-fields provided by ESO for 
each mask set-up. Large-scale intensity fluctuations were 
removed by dividing the flat-fields by a smoothed version 
using a 50 pixel median. Then the multi-object spectra 
were divided by the resulting normalised flat-fields. 

The wavelength calibration was carried out in two 
steps. First, the dispersion relation was derived using He, 
Hg/Cd and Ar spectra, which were obtained with the same 
mask configuration as the object spectra. These calibra- 
tion spectra were used to obtain a fourth-order dispersion 
relation for each row of the MOS frame. Typical relative 
uncertainties of the relation were ~ 0.5 A (~ 1/10 pixel). 

The calibration spectra were obtained in zenith po- 
sition during daytime. Therefore, bending and the posi- 
tioning of the laser-cut masks in the mask exchange unit 
(MXU) could cause variations of the slitlet locations be- 
tween calibration frames and object spectra of up to one 
pixel. This effect was corrected by measuring the position 
of the prominent telluric line [O I] A 5577 and by shifting 

1 Munich Image Data Analysis System 
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the spectra accordingly. The final error of the wavelength 
calibration was found to be about 0.5 to 1 A. 

The sky spectrum underlying the object spectrum was 
derived by interpolation and subtracted. The interpolation 
was carried out by polynomial fits of first to third order 
as well as median fits using 7 to 15 pixels. In each case 
the most suitable approach was chosen according to its 
success in reaching an undistorted object continuum com- 
patible with the photometric data and the minimisation 
of sky line residua. The uncertainty of the zero-flux level 
of the flux calibrated one-dimensional spectra was found 
to correspond to a few 10~ 23 Wm~ 2 A _1 in the optical 
wavelength region. 

For the extraction of onc-dimcnsional object spectra 
from the two-dimensional MOS spectra the S/N-optimised 
algorithm of Home ( 1986) was used. Pixels with a devi- 
ation greater than 5 a from the expected intensity profile 
were ignored in the variance-weighted summation proce- 
dure. This allows eliminating the majority of cosmics and 
CCD defects. As a by-product of the variance estimates 
the algorithm provides for each final spectrum /(A) an 
error function A/(A) representing essentially the photon 
noise level. 

The atmospheric extinction correction was carried out 
using the ESO standard extinction curve fTug ll9f7jl . 

A first-order flux calibration was performed using 
spectra of spectrophotometric standard stars (Oke 11997)1 
Hamuy et al. H992lll994f) usually obtained each night with 
5" wide slits. These spectra were reduced like the pro- 
gramme spectra and used to convert the ADUs into flux 
units. Then the flux-calibrated observed spectra were di- 
vided by the corresponding standard star spectra from the 
literature in order to obtain response curves. Eventually, 
representative correction curves were created for each run 
combining the individual response curves and smoothing 
the result by a spline interpolation. 

The limited slit width of 1" unavoidably causes flux 
losses in the case of extended objects, modest seeing 
and/or decentred slit positions. These flux losses were cor- 
rected using the total photometric fluxes measured during 
the photometric programme (Heidt el al. I2003b|) . In prac- 
tice, quasi-photometric fluxes were calculated from the 
spectra (applying the known FORS filter response curves) 
and divided by the photometric results to determine cor- 
rections. Usually the averaged factors of the g and R filters 
were used. For very blue or very red SEDs g and R were 
substituted or supplemented by B or /. 

The resulting spectroscopic fluxes can still be affected 
by inaccuracies of the atmospheric extinction law, the re- 
sponse curve, the zero-flux level and residua of strong sky 
lines. For extended objects with strong colour gradients 
the slit loss correction can be incorrect because of system- 
atic deviations between the measured and the represen- 
tative object spectrum. To estimate the error introduced 
by such effects, the spectroscopic and photometric fluxes 
of the final co-added spectra were compared. For the pri- 
mary spectroscopic targets this analysis indicates typical 



relative flux uncertainties of about 5 % in the g and R 
range and 15 % in the B and I range. 

The flux-calibrated spectra were S/N-optimised co- 
added using the squared ratios 



to 



(/(A)) 
(A/(A)} 



(1) 



of the suitably averaged spectra /(A) and error functions 
A/(A) as weights. At very low S/N errors of the spec- 
tral continuum level are no longer negligible, which may 
lead to unrealistic weights. Hence, an additional wo = 0.5 
was introduced in Eq. |JT]|. assuming that the noise levels 
of the individual spectra are of the same order. To re- 
move artifacts such as CCD defects or undetected cosmics 
a cr-clipping procedure was carried out before the averag- 
ing routine was started. The clipping threshold was set to 
6A/(A). 

Since the detector oversampled the spectra with about 
4 pixels per resolution element, all spectra were smoothed 
to match the spectral resolution using a Gaussian filter. 
The error functions were adjusted accordingly, converting 
noise per pixel into noise per resolution element. 

Finally, the spectra were corrected for Galactic dust 
extinction using the law of Cardelli et al. (1989) and the 
expected colour excess Eb-v — 0.018 mag of the extra- 
galactic objects in the FORS Deep Field (Schlegel et al. 
1998). In the visual the correction amounts to about 5 %. 

The atmospheric B (6830 - 6950 A) and A (7570 - 
7710 A) absorption bands were removed by dividing the 
spectra with a band model obtained from the average of 
a large number of high S/N spectra. 

The resulting data set of the FDF spectroscopic sur- 
vey contains 604 co-added, one-dimensional spectra, 51 % 
of which represent primary targets. Fig. [21 shows the S/N 
versus / magnitude for these spectra. Primary and sec- 
ondary objects are distinguished by different symbols (cir- 
cles and crosses). As expected, the average S/N is higher 
for the primary objects (« 20) than for the secondary ob- 
jects (w 4). 24 % of the spectra have a continuum S/N < 1. 
These objects, essentially all of which are secondary tar- 
gets, were excluded from the further analysis. 

5. Derivation of redshifts and object types 

Before analysing the spectra in detail, we derived for all 
objects redshifts and rough object types using the follow- 
ing iterative procedure: 

(1) An approximate redshift was determined from a visual 
line identification. 

(2) An approximate object type was determined by com- 
parison with known spectral energy distributions. 

(3) By averaging the spectra for each object type a set of 
empirical templates was constructed. 

(4) Improved redshifts and types were derived by cross- 
correlating the object spectra with these template 
spectra. 
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(5) Steps (3) and (4) were repeated until an optimal fit 
was achieved. 



In detail we proceeded as follows: 



5.1. Redshifts from visually identified lines 

For low-redshift galaxies with lines and blends with known 
rest wavelength, redshifts were determined (preferably us- 
ing strong emission lines) by fitting Gaussians to the line 
centres. Spectra of high-redshift galaxies are characterised 
by resonance lines and blends with complex profiles, such 
as blue-shifted absorption lines and P Cygni profiles (see 
e.g. Adelberger et al. I2003f) . Hence, the visual redshift 
derivation was based on features such as the emission lines 
He II A 1640 and C III] AA 1907, 1909 and the low ionisation 
absorption lines Sill A 1526, Fell A 1608. These lines were 
expected to represent the systemic velocities of the galax- 
ies or outflows not larger than a few lOOkm/s. Thus, our 
redshifts should be consistent with those of Pettini et al. 
(USUI, Frye et al. pUU2j) . and Adelberger et al. PUU3)1 . 



5.2. Spectral definition of object types 

Since most of the FDF galaxies are not resolved classifi- 
cation schemes based on the morphological Hubble types 
as applied by Kennicutt l|1992[) and Kinney et al. (1996) 
are not useful in our case. Therefore, we created a rough 
system of five galaxy spectral types (I to V) differing in 
the ratio of the rest-frame UV flux to the optical flux. The 
strengths and profiles of spectral lines have almost no im- 
pact on this classification. Obviously a one-dimensional 
sequence of only five types cannot be sensitive to the sub- 
tleties of the composition of the stellar population, the 
influence of dust and other characteristic galaxy parame- 
ters. However, in view of the greatly varying S/N of our 
spectra and the limited accessible wavelength range, our 
classification has the advantage of being rather robust to 
S/N induced effects. 

Our classification obviously requires the presence of 
reliable continua. Therefore, for the strong Lya emitters 
with weak continua a supplementary template VI was 
used. 

QSOs (or, more generally, galaxies dominated by an 
active galactic nucleus), whose spectra could not be fit- 
ted by the spectral standard types I to V were allotted a 
special class VII. 

Stars (type VIII) were identified on the basis of their 
angular intensity profiles and by comparison with stel- 
lar reference spectra from the literature (e.g. Pickles 
1985J. Their spectra were classified following the system of 
Morgan and Keenan (Morgan & Keenan fl973l Keenan & 
McNeil 1976) extended to very cool objects (e.g. Martin 
et al. lT555jl . 



5.3. Calculation of the template spectra 

The template spectra were calculated by averaging FDF 
spectra of the corresponding type with secure redshifts. 

As a preparatory step the individual spectra were 
masked to exclude wavelength regions affected by system- 
atic errors. The resulting useful wavelength range usually 
included the interval between 3800 and 8250 A, except 
for regions containing strong sky line residuals. Then the 
spectra were converted to the rest-frame system and re- 
binned with a step size of lA. 

To obtain average spectra with equivalent widths 
(EWs) equal to the average EW of the individual spectra, 
the masked rest-frame spectra J™ (A) were divided into a 
continuum and a line part according to 

/r(A) = / c m i (A)/ 1 -(A). (2) 

The f"\(X) were generated using a median filter of 100 A 
width. 

The /™(A) were averaged as follows: 



1 " 

where 

n 

M(A) =J2 "hW 

i=l 



(3) 



(4) 



is the sum of the mask functions m^(A) (values either or 
1) of all n spectra. 

The continua were added sorted by redshift in order to 
maximise the respective size of the overlapping wavelength 
range. Step-by-step the summation was realised using 



/"(A) 



(5) 



with F c (A) = 0, Cx = 1 and the effective flux correction 
factor for i > 2 



M i _ 1 (A)/ c m t (A) 

F c ,i-l(\) I fm Fc ._ 1 



1^0 



where 



Mi(A) = X>iW. 

3=1 



(6) 



(7) 



Ci ensures that all spectra are added with the same weight. 
Finally, the average continuum was calculated by 



1 



FJX) = — —F cn (\) 
V ' M(A) ' V ' 



(8) 



and suitably normalised. 

F C (A) was derived adding spectra in ascending as well 
as descending redshift order. The relative deviations were 
smaller than 5%. The final templates were calculated ac- 
cording to 



F(X) = \(F cA (X)+F cd (X))F 1 (\). 



(9) 
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Fig. 3. Empirical templates used for classifying the galaxies of the FDF spectroscopic sample. The noise level is 
indicated by a dashed line. SEDs I to V are normalised to the same mean flux in the wavelength range 3500 to 4500 A. 
SED VI for bright Lya emitters (panel (f)) is adjusted to have the mean flux of SED V between 1400 and 1600 A. For 
comparison SED V is overplotted in panel (f) (dotted line). 
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Fig. shows the resulting templates for the types de- 
fined in Sect. 15.21 As shown by the figure the different 
types differ mainly in their ratio between the UV and op- 
tical flux. However, the nebular emission line strengths 
of the templates also increase from SEDI (no emission) 
to SED V (numerous strong emission lines). Template VI, 
which was constructed from SED V spectra with strong 
Lya emission, shows a similar continuum as SED V. 

As they were composed from spectra of different red- 
shift the mean spectra cover large wavelength ranges. 
Since different spectral regions originate from objects of 
different redshifts, it is clear that the mean spectra may 
not correspond to real galaxy spectra, if evolutionary ef- 
fects are significant. Nevertheless, they are well suited for 
the purpose of deriving redshifts and rough object types. 

5.4. Redshift derivation 

To derive an accurate uniform set of redshifts, all spectra 
were (after conversion to a logarithmic wavelength scale) 
cross-correlated with their corresponding templates (see 
SimkinHnZH Tonry & Davis ll97(?| . Only regions of the in- 
put spectra with good quality were included in the calcula- 
tion. Furthermore, the continua of spectra and templates 
were removed to suppress the continuum background in 
the correlation function. 

The most likely redshift was derived using a \ 2 test 
checking the correspondence of the respective spectrum 
(continuum inclusive) and the template, with the redshift 
being a free parameter. This procedure takes into account 
that the amplitude of cross-correlation peaks depends on 
object type, redshift and S/N. 

5.5. Final redshifts and spectral types 

For each spectrum the redshift derivation procedure was 
carried out for each template shown in Fig. [3] Then the 
most probable redshift and spectral type were determined 
from the best fit. If the results for two templates indicated 
similar redshifts and fit quality, the average was taken. 
In general, the redshift differences were within the error 
limits (see below). 

Since low S/N or unusual SEDs can introduce errors 
in automatic procedures, all results of the automatic red- 
shift derivation were verified by visual inspection. In this 
way 341 redshifts could be confirmed to be obviously cor- 
rect (277) or very likely (about 90% confidence) correct 
(64). These objects are listed in Table They represent 
88 % of the primary and 23 % of the secondary targets. 
Considering the objects with secure redshifts only, the val- 
ues are 77% and 14%. For the remaining 263 objects the 
redshift derivation was either ambiguous or no redshift 
could be derived because of too low a S/N. 

To estimate the redshift accuracy we compared our 
results with those obtained from higher-resolution spec- 
tra available for part of our galaxies. For 51 brighter 
objects with z < 1, we made use of medium resolu- 
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Fig. 4. / magnitudes as a function of the spectroscopic 
redshifts of the extragalactic Catalogue objects. Primary 
and secondary targets are indicated by circles and crosses, 
respectively. In contrast to Fig.^ the special primary ob- 
jects discussed in Sect. [5] are marked by triangles. The 
bright object with z = 3.365 is the quasar Q 0130-260. 



tion FORS-600R spectra (AA ~ 4 A) obtained for a dif- 
ferent programme (see Ziegler et al. 2002 Bohm et al. 
2004J. The comparison indicates relative redshift errors 
Az/(z + 1) ~ 4 x 10~ 4 for the 150 1 spectra showing emis- 
sion lines. For type I spectra, which are characterised by 
relatively weak absorption features, the accuracy is about 
8 x 10~ 4 . 

To estimate the redshift uncertainties for all spectra, 
especially the high-redshift ones, the dependence of the 
redshift on the template used was taken into account. 
Moreover, redshifts were also derived by measuring spec- 
tral line positions. Taking together the different error es- 
timates the typical relative uncertainties range between 
4 x 10 -4 and 10 x 10~ 4 depending on redshift regime, 
spectral type and S/N. Values for the individual redshift 
errors are presented in Tabled 



6. The Catalogue 

The basic properties of all 341 objects with certain or 
probable redshift are listed in Table (in the following 
referred to as 'the Catalogue'), which is available in elec- 
tronic form only, except for a sample page presented in 
Table The individual columns of the Catalogue have 
the following content: 



No.: Object number according to Heidt et al. I)2003bj) . For 
a few objects not listed in the photometric catalogue 
(mainly because no photometry could be obtained due 
to crowding) new numbers, starting with 9001, have 
been assigned. 

RA: Right ascension (equinox: 2000) in hours, minutes 
and seconds. 

DEC: Declination (equinox: 2000) in degrees, arcminutes 
and arcseconds. 
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Table 1. Excerpt from the electronically available Catalogue of the FDF spectroscopic sample. Descriptions of the 
individual columns are given in the text. 



No. 


RA 


DEC 


B 


/ 


T 

' cxp 


f/fo 


S/N 


Q s 


Type 


z 


dz 


a 


Q z 


Rem. 




01:... 


-25:. . . 


[mag] 


[mag] 


[min] 


[%] 


(Band) 










[%] 






4522 


06-03 8 


45:42 


99.0 


25.4 


1121 


97 


3 4 (I) 


2 


v 


4.9996 


0.0048 


43 


2 




4573 


06:04.0 


46:55 


25.5 


23.9 


281 


31 


3.6 (R) 




IV 


6993 


0.0007 


41 


2 




4609 


06:04.2 


46:45 


25.3 


23.4 


135 


47 


4 1 (a) 




IV 


0.6431 


0.0007 


50 


1 




4654 


06:04.3 


47:19 


22.2 


20.1 


30 


5 


4.0 (R) 




IV 


3989 


0006 


40 


2 




4657 


06:04.3 


43:20 


21.8 


19.3 


120 


2 


8 8 (a) 

o.o \y ) 


1 


III 


0.2244 


0005 


36 


1 


600R 


4662 


06:04.3 


43:12 


25.4 


24.0 


240 


35 


3 5 (a) 

o.o yyj 


2 


v 


0.7386 


0.0007 


50 


2 




4667 


06:04.3 


47:14 


22.6 


20.7 


30 


16 


7.5 (i?) 




IV 


0.4542 


0006 


32 


1 




4682 


06:04.4 


46:15 


99.0 


23.6 


510 


49 


6 1 (a) 


2 


v 


3.1428 


0033 


18 


1 




4683 


06:04.4 


46:51 


20.2 


18.6 


155 


54 


284 2 (a) 




VII 


3 3650 


0.0070 




1 


QSO 


4684 


06:04.4 


48:00 


25.7 


23.6 


30 


87 


2 7 (a) 
z " 1 \y } 




IV 


0.7843 


0.0007 


52 


2 




4685 


06:04.4 


43:11 


25.0 


23.5 


240 


50 


9 4 (a\ 

\y> 




IV 


2.6254 


0025 


16 


1 




4686 


06:04.4 


45:54 


22.7 


21.8 


442 


54 


43 8 (n) 
<±o.o yyj 




IV 


0.0943 


0.0004 


15 


1 




4691 


06:04.4 


48:34 


25.8 


24.3 


550 


66 


13 1 (a) 
±o.± \y ) 




V/VI 


3 3036 


0.0043 


67 


1 


LAB 


4701 


06:04.5 


42:52 


29.4 


24.7 


225 


60 


1 4 (T) 




v 


4.9132 


0.0047 


82 


2 




4729 


06:04.6 


47:36 


22.3 


21.1 


120 


51 


34 2 (a\ 
o'i-^ \yj 




VIII 


o 






1 


G 


4745 


06:04.6 


44:03 


25.2 


23.8 


434 


67 


10 2 (a) 
iu,i \y i 




v 


2.6173 


0029 


17 


1 




4752 


06:04.7 


46:54 


26.5 


99.0 


473 


29 


3-2 (g) 


2 


V/VI 


3.3802 


0.0044 


119 


1 


LAB 


4795 


06:04.8 


47:14 


24.4 


23.3 


570 


55 


15.9 (g) 




IV 


2.1593 


0.0022 


14 


1 




4871 


06:05.1 


46:04 


24.9 


23.4 


490 


58 


15.4 (g) 




IV 


2.4724 


0.0024 


13 


1 




4882 


06:05.1 


48:38 


21.3 


19.0 


120 


31 


75.0 (R) 




II 


0.2777 


0.0006 


11 


1 




4910 


06:05.2 


46:04 


25.5 


24.3 


120 


40 


2-1 (S) 




V 


0.8437 


0.0007 


65 


2 




4954 


06:05.4 


46:20 


25.8 


23.8 


160 


37 


2.1 (R) 




IV 


0.9703 


0.0008 


55 


2 




4993 


06:05.5 


45:56 


26.5 


21.8 


245 


50 


11.7 (/) 




I 


0.7632 


0.0014 


24 


1 




4996 


06:05.5 


46:28 


24.4 


23.3 


281 


56 


14.6 (g) 




IV 


2.0278 


0.0021 


15 


1 





B: Total apparent B magnitude (Vega system, see Hcidt 
et al. I2fl03b| . Non-detections in B are marked by the 
value 99.0. 

I: Total apparent I magnitude (as above). 

T CX p: Total exposure time in minutes. 

f / fo'. Ratio between the flux which passed through the 
slit and the actual object flux in %. Low f/fo usually 
correspond to large object extensions. Typical values 
for point-like objects are around 70%. Large values 
(> 80%) can be caused by systematic spectral errors 
and/or inaccurate photometry due to very low fluxes 
or object crowding. 

S/N: Average signal-to-noise ratio per resolution element 
in the filter band given in parentheses (B, g, Roy I). In 
each case the band with the highest S/N was selected. 
The S/N as a function of wavelength was calculated 
by dividing the object spectrum by its error function. 

Q s : Flag indicating systematic errors in the spectrum. 
Q s — 1 refers to distorted spectra, Q s = 2 to local 
defects. 

Type: Object type as defined in Sect. 15.21 and Fig. [31 The 
classes I to VI represent galaxies, VII indicates QSOs, 
and VIII refers to stars. 

z: Spectroscopic redshift. 

dz: Mean error of the redshift. 

a: Relative rms deviation between spectrum and the op- 
timal template in % of the average spectral flux. 



Q z : Quality of the redshift. Q z = 1 indicates objects with 
secure redshifts and Q z = 2 with probable redshifts 
(90% confidence level). 

Rem.: Further information on the object. For stars a 
rough spectral type is given. Quasars and strong Lya 
emission galaxies are indicated by the entries 'QSO' 
and 'LAB' (Lya bright, i.e. Lya emission EW > 20 A), 
respectively. '600R' indicates galaxies whose redshift 
and object type were verified by means of the spec- 
troscopic data of the medium resolution spectra (see 
Sect. 151)1) . 

To illustrate the sample properties Fig. shows the 
I magnitudes of the extragalactic Catalogue objects as a 
function of spectroscopic redshift. Primary (236) and sec- 
ondary objects (63) are indicated. High-redshift primary 
targets, which were selected in a special way (12 objects, 
see Sect.EJ), are marked by triangles to distinguish them 
from the I magnitude selected primary targets, which are 
indicated by circles. 



7. Basic properties of the spectra 

In addition to the catalogue the flux-calibrated individual 
object spectra with reliable redshifts and no strong sys- 
tematic errors are also available as an electronic atlas. In 
the following we describe some properties of these spectra. 
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Fig. 5. Examples of the FDF spectra. The positions of some prominent spectral lines are marked by vertical dotted 
lines. In panel (b) the position of the UV absorption feature at 2200 A is indicated by a horizontal bar. Vertical 
dash-dotted lines indicate the positions where terrestrial atmospheric features ([OI] A 5577, NalD, [OI] AA 6300, 6364, 
B band, and A band) were removed. The noise level is indicated by a dashed line. At the top of each diagram the 
object number, the spectral type, and the rcdshift are given. 



7.1. Stars 

The FDF spectroscopic sample contains 42 stars of spec- 
tral types G to L. Most frequent are K stars (18). Five 
stars have spectral types later than M5. None of the stars 
is brighter than 17 mag in /. Since the FDF is located 
close to the South Galactic Pole, most of the observed 
stars must belong to the halo or the thick disk. 



7.2. Galaxies 

The spectra of the 158 galaxies with z < 1 show a great va- 
riety of spectral types ranging from ellipticals to extreme 
starburst galaxies. 

Since early-type galaxies are intrinsically red, it is dif- 
ficult to find such objects with high redshifts in optical 
surveys. For z > 1.3 the 4000 A break is shifted beyond 
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our spectral window. This readily explains why no type I 
galaxy beyond z = 1.17 (FDF-2256, see Fig. EJa)) could 
reliably be identified in the FDF, so far, and why only 
six candidates of spectral type I or II with z > 1 were 
detected. 

Some of the intermediate redshift galaxies show a de- 
pression at about 2200 A (rest frame). As an example 
FDF-4049 is presented in Fig.GJb). The 2200 A feature of 
this galaxy corresponds to Sb-v ~ 0.15 mag if a Galactic 
extinction law is assumed. The figure seems to illustrate 
that the well-known broad dust absorption feature, which 
is thought to be caused by small carbon-rich particles (see 
Cardelli et al. H989[) . is also present in at least part of the 
intermediate-redshift galaxies. 

In the rest-frame wavelength range between 2000 and 
3700 A galaxy spectra are poor in prominent spectral lines. 
This complicates the derivation of reliable redshifts in the 
range between 1.5 and 2. Therefore, in this redshift range 
secure redshifts could be derived for two objects only. 

For z > 2 we observe only starburst galaxies. Their 
spectra are dominated by the Lya line (with profiles rang- 
ing from strong absorption to strong emission), the ab- 
sorptions of the Lya forest short-ward of 1216 A, and the 
spectral break of the Lyman limit at 912 A. Furthermore, 
the spectral region red-ward from Lya is marked by nu- 
merous lines, usually showing absorption profiles, and of- 
ten of interstellar origin (e.g. lines of Si II and C II). Except 
for Lya significant emission could be detected in a few 
cases only, either in the form of a P Cygni profile (e.g. 
C IV AA 1548, 1550) or as pure emission (e.g. He II A 1640 
and C III] AA 1907, 1909). All these features are visible 
in the spectra of FDF-5903 (Fig. [He)) and FDF-6024 
(Fig.Etd)). 

Four (44%) of the nine galaxies with z > 4.5 were 
found to be Lya bright (type VI). This percentage is com- 
parable to the redshift range 3 < z < 4 where about 37 % 
are LABs. Up to now, the objects with the highest spec- 
troscopic redshifts in the FDF are the galaxies FDF-4522 
and FDF-5812 (Fig.E^e)) both having z « 5.0. 



7.3. Quasars 

Most of the eight quasars at redshifts between 0.865 and 
3.365 identified in the FDF show normal QSO spectra 
dominated by broad emission lines of Lya, NV, Si IV, 
C IV, C III] and/or Mg II. The only exception is FDF-6233 
at z = 2.32 (see Fig.GJf)) which belongs to the rare class 
of broad absorption line quasars (BAL, e.g. Menou et al. 
2001). Its spectrum is dominated by strong, broad, blue- 



shifted absorption troughs of Si IV, C IV and Al III indi- 
cating outflows of several 1000 km/s. Lya, N V, and CIII] 
show relatively narrow emission components. In contrast 
to the other quasars the continuum of FDF-6233 decreases 
rapidly short-ward of 1750 A (rest frame). 
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Fig. 6. Redshift distribution of the FDF spectroscopic 
sample for (a) the complete Catalogue (see Sect. EJ), (b) 
- (f ) galaxies of different object type (see Fig. |3J) and (g) 
the quasars. In diagram (f ) galaxies of type V (white bars) 
as well as type VI (black bars) are plotted. The redshift 
resolution is Az = 0.05. 



8. Distribution of spectroscopic redshifts 

Fig. Efa) shows the spectroscopic-redshift distribution of 
the galaxies and QSOs with a resolution of Az = 0.05. 
Significant local maxima are evident at z = 0.325 (29 ob- 
jects), z = 0.775 (18), z = 2.375 (11) and z = 3.375 (10). 
In part this clustering of redshifts had already been no- 
ticed in the distribution of the photometric redshifts. In 
some cases the overdensities extend over more than two 
bins. For instance, the three adjacent bins comprising the 
redshift range between 2.30 and 2.45 contain 26 objects 
which represent about 27 % of the z > 2 sample. On the 
other hand, the peak around z = 0.325 is mainly pro- 
duced by a redshift interval ranging from 0.331 to 0.339 
((z) — 0.335, <t z — 0.002) which is smaller than the bin size 
but which includes 23 galaxies. As indicated in Fig.[7]such 
fine structure of the redshift distribution is also present at 
high redshifts. An example is the narrow redshift range 
between 2.372 and 2.378 ((z) = 2.374, <j z = 0.002) which 
includes eight galaxies, or about 73 % of the objects within 
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Fig. 9. Distribution of galaxies within the range 3.36 < 
z < 3.41 in the FDF. The position of the quasar Q 0130- 
Fig. 7. Redshift distribution of the FDF spectroscopic 260 is indicated by a circle. The three objects close to the 
sample in the range 2 < z < 3.5 plotted with a redshift quasar are strong Lya emitters, which were selected using 



resolution of Az = 0.01. 



narrow-band images (see Sect. • 



-25.72 


1 


1 i 1 i 
i' + 


i 1 i 1 i 1 i 1 i 
0.33 < z < 0.34 " 


-25.74 




+ 

+ 




-25.76 






+ + 


-25.78 




+ 


+ 

+ 


-25.8 




+ 

4- 


+ + 


-25.82 


r 


i 1 i 1 


+ 





16.58 16.56 16.54 16.52 16.5 16.48 16.46 
RA [deg] 



Fig. 8. Distribution of galaxies within the range 0.33 < 
z < 0.34 in the FDF. 

the standard bin around z = 2.375 (see above). As pointed 
out by Frank et al. I|2003p the high-redshift galaxy over- 
densities correlate well with QSO metal absorption line 
systems in the direction of the FDF. Besides the striking 
peaks, there are also indications of marked gaps in the 
redshift distribution. At low redshifts this is particularly 
evident around z = 0.525 where the redshift distribution 
reaches a local minimum. At high redshifts a lack of galax- 
ies exists between 2.80 and 3.05, where no objects could 
be detected, so far. The relatively low number of objects 



between 1.5 and 2.0 is a consequence of the difficulty to 
identify galaxies in this redshift range from visual spectra 
(see Sect. EH - 

Since the spectroscopic sample is representative for the 
most luminous objects at each redshift (see Sect. El, the 
observed distribution can (apart from the redshift range 
1.5 < z < 2.0) be assumed to be representative of the 
real distribution for the absolutely bright objects in the 
direction of the FDF. On the other hand, we have almost 
no information on the distribution of the intrinsically faint 
z > 2 objects (Mb > —20 mag), as such objects are rare 
in our sample. 

Since the first photometric observations of the FDF a 
(poor) cluster is known to be located in the southwest- 
ern corner of the field (see Heidt et al. I2003bl Ziegler 
et al. I2004J) . This cluster contributes to the peak of the 
FDF redshift distribution at z w 0.335. As indicated by 
Fig- |H1 the other objects at this redshift are not randomly 
distributed cither, as no single object of this redshift is 
found in the northwestern corner of the FDF. For the 
other peaks no significant deviation from a random dis- 
tribution over the area of the FDF could be detected. 
This negative result may in part be due to the insuffi- 
cient number of objects. But, since the angular size of 
the FDF covers just about the size of a large galaxy clus- 
ter at high redshift no major angular clustering effects 
were expected. Nevertheless, some indications for cluster- 
ing on small scales were found. For instance, the galaxies 
between z = 3.36 and 3.41 populate the southern part 
of the FDF only (Fig. |5J). Furthermore, three Lya emit- 
ters of this redshift were found less than 5" from the 
bright quasar Q 0130-260. Moreover, several close pairs 
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Table 2. Distribution of objects with reliable redshifts 
in the FDF spectroscopic sample. The table indicates for 
each type the total number of objects N, the median red- 
shift z, the average redshift (z) and its variance a z . 



Type 


N 


z 




o z 




I 


32 


0.4 


0.6 


0.3 


°< 


II 


25 


0.3 


0.5 


0.3 


IS 


III 


50 


0.6 


0.7 


0.4 


o 


IV 


96 


1.0 


1.4 


1.0 


o 

lO 


V + VI 


88 


2.4 


2.4 


1.3 




QSOs 


8 


2.3 


2.3 


0.8 












o 


Extragalactic 


299 


0.9 


1.4 


1.2 





objects 
Stars 



42 



of galaxies with similar redshifts were discovered. In one 
case a dense group at z « 2.347, comprising four members 
(FDF-5135, 5165, 5167, 5190), was identified. The mutual 
distances of these objects are < 3". 

In general, the FDF spectroscopic data, indicating pre- 
dominantly moderate clustering of galaxies up to high red- 
shifts, which can be attributed to large-scale structure, 
agree well with earlier results of, e.g., Adelberger et al. 
(UnnSJ, Steidel et al. fiMSj l. or Cohen et al. lj2TOjl . 

The redshift distribution of the different object types 
defined in Sect . 15 . 21 and 15 . 31 is shown in Fig.HJb)-(g). The 
type-dependent median and mean redshifts are presented 
in Table FJ1 For the reasons outlined above the types I 
and II are concentrated at low redshifts (median values 
0.4 and 0.3). Furthermore, type I shows a strong over- 
density at z ~ 0.335 where 11 of 32 early-type galaxies 
are located. This obviously reflects the predominance of E 
and SO galaxies in the cluster at this redshift. The change 
of the apparent composition of object types with redshift 
due to the redshift-dependent selection bias is particularly 
obvious for class V/VI: At low redshifts this type plays a 
minor role (7% of the observed objects at z < 0.5), but at 
high redshifts type V/VI becomes the dominating class. 
At 2 < z < 3 the extreme starburst galaxies contribute 
38 % to all identified objects, at 3 < z < 4 they present 
78 % and at z > 4 the fraction increases to 91% (10 of 11 
objects). 



9. Exploring evolutionary effects 

A detailed discussion of the whole FDF spectroscopic sam- 
ple is beyond the scope of this paper. Instead, we will focus 
on some properties of the high-redshift objects. As aver- 
age spectra of spectroscopic subsamples selected by red- 
shift and/or spectral type are known to be well suited to 
investigate the characteristic properties of galaxy popula- 
tions, the following discussion is based mainly upon such 
mean spectra. An analysis of individual object spectra can 
be found in, e.g., Mehlert et al. l|2002fl and Tapken et al. 
(2001. 
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Fig. 11. Luminosities at 1500A of the z > 1.35 FDF 
galaxies as a function of redshift. To derive the luminosity 
the spectral flux was measured in the wavelength interval 
1480 to 1520 A. Galaxies used for the calculation of the 
composite spectra are indicated by circles. The remaining 
high-redshift objects are marked by crosses apart from 
QSOs, which are indicated by asterisks. 



9.1. The calculation of composite spectra 

Composite spectra of FDF subsamples were calculated us- 
ing similar procedures as used bv lS03l to facilitate a com- 
parison with their spectroscopic sample (see Sect. l9.4|) . In 
detail, the flux-calibrated, co-added spectra (see Sect, rj} 
were shifted to the vacuum rest frame, rebinned to a scale 
of 1 A per pixel, scaled to a common mode in the wave- 
length range between 1250 and 1500 A, and finally aver- 
aged unweighted to get a representative average. 

The averaging of spectra showing different continuum 
levels may lead to a systematic change of the mean equiva- 
lent width of spectral lines. A comparison with mean spec- 
tra based on the alternative method described in Sect. 15.31 
showed that for Lya an overestimate of the true EW 
of the order 10% can occur. For other lines no signifi- 
cant effect could be detected. The dependence of the Lya 
strength on the normalisation method must be kept in 
mind when analysing the results described in the follow- 
ing (see Sect. I9.5|) . 



In contrast to IS03I we did not exclude extreme flux 
values during the averaging. Instead, our averaging pro- 
cedure omitted wavelength ranges affected by significant 
sky line residuals. 
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Fig. 10. Composite spectrum of 64 FDF galaxy spectra with S/N > 4 in the range 2 < z < 4 (thick solid line). The 
positions of prominent UV lines are indicated. The thin solid line at the bottom of the spectrum gives the mean error 
of the average spectrum as derived from the scatter of the individual spectra. Outside strong spectral features this 
error agrees with the value expected from the observational errors. In the case of some strong features such as Lya 
the mean error is dominated by the physical variance in the individual spectra. 



9.2. The overall mean spectrum of the high-redshift 
sample 

In Fig. El we present a composite of 64 galaxy spectra 
with 2 < z < 4 and a continuum S/N > 4. This sub- 
sample of photometrically selected objects (I < 24.5 mag) 
represents about 91 % of all galaxies with secure redshift 
in the range 2 < z < 4. All objects are primary targets 
(see Sect. 01. Three z > 3 galaxies (FDF-7452, 7644, 8215) 
were initially included in the observations because of sus- 
pected Lya emission. One z > 3 object (FDF-4691) was 
originally a secondary target showing strong Lya emis- 
sion. One z < 3 object (FDF-9011) was included due to 
its location close to another sample galaxy. The individ- 
ual spectra are representative of the integrated light of 
the entire galaxies, since the slit losses are not more than 
twice those for point-like objects (see Sect.EJ. The mean 
redshift of this subsample is (z) = 2.72 (cr z = 0.49). The 
UV luminosities at 1500 A of the selected galaxies (derived 
from the flux in the interval 1480 to 1520 A) are indicated 
by circles in Fig. El (log Li 50 o[W/A]} amounts to 34.22 
(ff£ = 0.22). 



The composite spectrum reaches a continuum S/N rs 
60, which allows a comparison with similar high-S/N high- 
redshift spectra from the literature (e.g. Pettini et al. 2000 
S03). While the general character of all these spectra is 
similar, differences in spectral details, such as the Lya 
strength and the continuum slope, are evident. We discuss 
these differences in Sect. 19.41 

9.3. The redshift dependence of spectral properties 

To investigate evolutionary effects in our sample we com- 
puted in the next step mean spectra of the galaxies within 
the redshift intervals 2 < z < 3 and 3 < z < 4. The re- 
sulting composite spectra with mean redshifts (z) = 2.39 
(<j z = 0.19) and (z) — 3.33 (er 2 = 0.20) are based on 
42 and 22 objects, respectively. The results are given in 
Fig. Eland Table El 

Following the normal spectroscopic convention in 
Table (and throughout this paper) for net absorption 
lines we give positive equivalent width values, while net 
emission features are characterised by negative EW fea- 
tures. The EW mean errors were estimated on the ba- 
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sis of the mean errors of the composite spectra, which 
were derived from the scatter of the continuum-subtracted 
individual spectra, and the uncertainties concerning the 
continuum level determination. In contrast, the mean er- 
rors of the continuum slope j3 (calculated in the range 
1200 to 1800 A assuming /(A) oc \@ following Leitherer 
et al. l|2002|l ) were derived from the scatter of the values 
measured in the individual galaxy spectra (see Fig. I21|) . 
Usually the j3 of the mean spectrum and the average j3 
of the individual spectra were roughly the same, differing 
by less than 0.1, which shows that the composite spectra 
are characteristic for the respective samples, at least with 
regard to the continuum slope. 



As shown by Fig. El and Table the 2 < z < 3 mean 
spectrum has a rather flat continuum (/? = —0.56 ± 0.11) 
and Lya emission is weak (Wh ya = +7 ± 2 A). In con- 
trast, the 3 < z < 4 mean spectrum shows a steep UV 
continuum (/3 = —1.79 ± 0.13) and strong Lya emission 
(Why a = —10 ± 6 A). The equivalent widths of the ma- 
jority of spectral lines and in particular all interstellar ab- 
sorption lines show minor differences only. However, in 
general, the EWs tend to be higher in the lower redshift 
range. A strong redshift dependence is observed for the 
CIV resonance doublet, which shows a weaker absorp- 
tion at higher redshifts (as noted already by Mehlert et 
al. 12002). Furthermore, C III] displays a stronger emis- 
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Table 3. Properties of the FDF composite spectra repre- 
sentative of the redshift ranges 2 < z < 3 and 3 < z < 4 
in comparison. The table lists the number of spectra aver- 
aged, the mean redshift, the average luminosity at 1500 A 
(see Sect. I9.2J) . the Lya equivalent width, the continuum 
slope P, the mean equivalent width of low ionisation in- 
terstellar absorption lines (see Fig. I16f) and the EWs of 
CIV AA 1548, 1550, He II A 1640 and C III] AA 1907, 1909. 
All errors given are mean errors (see Sect. I9~T|) . 



z range 
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iogi 1500 [w/A] 


34.14 ±0.03 


34.35 ± 0.04 
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sion in the 3 < z < 4 mean spectrum. However, due to 
the line's location in the region of strong OH bands at 
z > 3 the 3 < z < 4 spectrum contains only five spectra 
at CIII]. Therefore, the observed increase of CIII] with 
redshift may not be statistically meaningful. 

The composite spectrum of the higher redshift range 
(3 < z < 4) includes only few objects with luminosities 
lower than logLi 50 o = 34.2 (3 out of 22; see Fig. ITl"|) . 
whereas 28 out of 42 objects in the range 2 < z < 3 
are below this value. To check the effect of the luminosity 
difference on the spectral properties, we also calculated 
composite spectra with logLisoo = 34.2 as lower lumi- 
nosity limit. The results showed that the luminosity does 
not affect the spectral differences, except perhaps for a 
slight decrease of the Lya strength in the high-luminosity 
subsample (from — 10 ± 7 A to — 5 ± 6 A), caused by the 
omission of two strong Lya emitters. 

9.4. Comparison with the composite spectrum of 
Shapley et al. 

S03 have published a composite spectrum of a homoge- 
neous sample of almost 1000 Lyman-break galaxies at 
z ~ 3 selected via the Lyman-break technique (see Steidel 
et al. l2003[) . Because of the different selection criteria, the 
redshift distributions of the FDF galaxies ((z) = 2.72, 
a z = 0.49) and the S03 sample ((z) = 2.96, a z = 0.29) dif- 
fer. 53% of the FDF high-redshift galaxies are at z < 2.5, 
whereas this redshift range is essentially absent in the 
S03 sample. Hence, in view of the results of Sect. 19.31 
differences in the spectral properties are to be expected. 
Therefore we calculated a composite spectrum of 28 FDF 
galaxies between z — 2.5 and 3.5 with a similar red- 
shift coverage {(z) — 3.11, a z = 0.28) as the S03 sam- 
ple. Fig. 1131 compares this composite spectrum with the 
S03 spectrum. The Lya strength ( Wh ya = — 9 ± 5 A) as 



Fig. 14. Redshift distribution of FDF galaxies with dif- 
ferent Lya types (see legend) . The presence of a quasar at 
the corresponding redshift is marked by an asterisk. 

well as the continuum slope ((3 = —1.54 ± 0.13) of the 
2.5 < z < 3.5 spectrum show good agreement with the S03 
values (Wi^ya = —11 A and (3 = —1.81). This shows that 
samples with comparable redshift distributions have sim- 
ilar average spectra independent of the selection method. 

The fact that the average apparent brightness of the 
S03 sample {R\ cga « 24.38 mag, converted from -Rab) is 
somewhat lower than that of the 2.5 < z < 3.5 FDF 
sample (i?vcga = 24.19±0.11mag) apparently has no effect 
on the average spectra. This is confirmed by a comparison 
of the composite spectra of the bright ({R) — 23.80 mag) 
and the faint half ((R) = 24.58 mag) of the 2.5 < z < 3.5 
sample, which show very similar spectra. The only possible 
difference may be a moderately higher Lya emission of 
— 16 ± 10 A for the faint subsample in comparison to —3 ± 
4 A for the bright subsample. 

9.5. Relations between Lya emission and other 
spectral properties 

In part of our spectra the Lya line occurs in emission. 
The emission originates from recombination of hydrogen 
ionised either by hot stars or by an AGN. The escape of 
Lya photons is hampered by multiple resonance scatter- 
ing resulting in large path lengths through the galaxy (see 
Chariot & Fall 1993). Hence, the Lya emission is very sen- 
sitive to the dust content and the geometry and structure 
of the HI velocity field (see e.g. Kunth et al.|:1998). 

To study the relation between the Lya emission and 
other spectral features we divided the high-redshift galaxy 
sample spectra into four classes according to the EW of the 
Lya emission component W^ ya ^. The four classes (T^ ya ) 
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Fig. 13. Comparison of composite spectra representative of the FORS Deep Field in the redshift range 2.5 < z < 3.5 
(thick line) and the S03 sample at z ~ 3 described in Steidel et al. (20021) and lS03l fthin line). The original S03 mean 
spectrum was smoothed to the same spectral resolution and normalised to the same mean continuum level in the 
interval 1250 — 1500 A. The lower panel gives the ratio between the FDF and the S03 composite spectra. 



were defined in the following way: 
Thy a = 0: No (detectable) emission. 
Thya — 1: Very weak emission (not measurable). 
T hya = 2: 2 A < -T4Wo < 20 A. 



Lya,c 



r Lya = 3: -^Lya.o > 20A. 

Mean and median rcdshifts of the individual Lya 
classes are listed in Table 0] The resulting redshift dis- 
tribution of the different types is shown in Fig. ^] The 
histogram indicates a conspicuous evolution of the typ- 
ical Lya characteristics from z — 2 to 5. At redshifts 
2 < z < 3 only 26 % of the objects shows significant Lya 
emission (TLya > 2). Only 3 of 18 Lya-bright galaxies 



(LABs, Thy a = 3) are at z < 3. They represent 6% of all 
galaxies in the range 2 < z < 3. At higher redshifts the 
fraction of galaxies with TLya > 2 increases to 67 % at red- 
shifts 3 < z < 4 and 73% at 4 < z < 5. Considering the 
galaxies with strong Lya emission only, the corresponding 
values are 37% and 36%. Hence, the evolution of the av- 
erage Lya strength between z = 3.3 and 2.4 is caused by a 
change of the frequency and the strength of Lya emission, 
if present. 

For the analysis of the correlation of Lya strength with 
other spectral features, we again took the 64 galaxy spec- 
tra in the range 2 < z < 4 and calculated composite spec- 
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Fig. 15. Comparison of FDF composite spectra representative of different Lya types ranging from strong emission 
(TLya = 3, uppermost spectrum) to pure absorption (TLya = 0, lowest spectrum). The spectra are normalised and 
shifted by multiples of 1.0 relative flux units. 



tra for each of the four Lya classes introduced above. In 
Fig-EDthe resulting spectra are plotted. To investigate the 
redshift dependence of the correlations, four further com- 
posite spectra were calculated, differing in redshift range 
(2 < z < 3 and 3 < z < 4) and Lya class T^ ya (0 — 1 and 
2 — 3) . To avoid too small numbers of spectra in individual 
composite spectra, the Lya types were grouped into two 
categories only. For all composite spectra EWs and contin- 
uum slopes were measured, allowing us to check whether 
the FDF sample follows the relations found by S03 



As pointed out by Pettini et al. ( 2000) and S03 promi- 
nent low-ionisation absorption lines of interstellar origin in 
spectra of high-redshift galaxies (like OI/SiIIA1303 and 



Si II A 1526) are normally saturated. The EWs of such lines 
are mainly determined by the product of the covering frac- 
tion of the neutral gas clouds and the velocity field, which 
makes them suitable indicators of the geometry and kine- 
matics of neutral gas in high-redshift galaxies. Fig. ^| 
shows the average EW of six prominent low-ionisation fea- 
tures as a function of the Lya equivalent width W^ ya . 
The four data points of the different Lya types (open cir- 
cles) indicate a significant decrease of Wlis with increasing 
W^Lya ■ The spectrum representing the strong Lya emitters 
shows a Wlis more than three times weaker than does 
the spectrum with pure Lya absorption, in close agree- 
ment with the results of Shapley et al. (crosses). Moreover, 
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Table 4. Distribution of Lya types T^ ya (as defined in 
Sect. 19.51) among the FDF high-redshift galaxies (2 < z < 
5). For each Lya class the total number of objects N, 
the median redshift z, and the mean redshift (z) with its 
standard deviation a z are given. 
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Fig. 16. The dependence of the strength of low-ionisation 
lines on the Lya equivalent width. W^lis is the average 
EW of the six strong low-ionisation features Si II A 1260, 
OI/SiIIA1303, C II A 1334, Si II A 1526, Fell A 1608, and 
Al II A 1670. The individual symbols represent the values 
measured in the FDF composite spectra for different Lya 
strengths and redshift ranges. Open circles show the values 
for the redshift range 2 < z < 4. Filled squares and trian- 
gles represent the smaller ranges 2 < z < 3 and 3 < z < 4, 
respectively. The 1 a error bars include the variance be- 
tween the individual spectra and the errors of the con- 
tinuum level definition. For comparison the corresponding 
data from IS031 are also included (crosses). 



significant differences of the relations for 2 < z < 3 
(filled squares) and 3 < z < 4 (filled triangles) cannot 
be detected, which suggests that a redshift dependence, 
if present, must be weak. Consequently, Fig. [TBI suggests 
a strong dependence of the Lya strength on the geome- 
try and kinematics of the H I clouds (in accordance with 
Kunth et al. [1998) independent of the redshift. 

Fig. IT71 shows CIVAA1548, 1550 as a function of the 
Lya strength. As CIV originates primarily in photo- 




Fig. 17. The dependence of the strength of CIV on the 
Lya equivalent width. The smaller symbols are as in 
Fig. ^3 The big filled square, filled triangle, and cross 
show the EWs measured in the 2 < z < 3, 3 < z < 4, and 
the S03 composite spectrum, respectively. 

spheres and winds of luminous hot stars (Walborn et al. 
1995; Heckman et al. I1998J) and the wind power depends 
on the metallicity, Wc iv is a good metallicity indicator 
(Leitherer et al. l2001f) . The figure confirms the strong de- 
pendence of the CIV strength on the redshift, noted al- 
ready by Mehlert et al. H2002f) . On the other hand, only 
a weak dependence of Wciv on WLya is indicated. The 
large Wciv mean errors suggest a marked scatter of the 
C IV strength for each Lya type. Consequently, the Lya 
emission appears to be rather independent of metallicity- 
dependent stellar wind properties. 

In Fig. ^] we investigate the dependence of the con- 
tinuum slope j3 on W\, Ya . As discussed in Sect. 19. 6| es- 
sentially traces the attenuation of hot star continua by 
interstellar dust. The diagram reveals a conspicuous red- 
shift dependence of the (3 — Wh ya relation. The galaxies 
at z < 3 (filled squares), which represent the main contri- 
bution to the FDF 2 < z < 4 sample (open circles), show 
a strong dependence of the continuum slope and Lya. In 
contrast, the galaxies at z > 3 (filled triangles) indicate 
a rather weak — Wh ya correlation, resembling the rela- 
tion for the S03 sample (crosses). The S03 (3 values had 
to be corrected (see Fig. I18|) to allow a comparison with 
the FDF values. Therefore, we computed Eb-v following 
S03 (using Starburst99 of Leitherer et al. (1999) and as- 
suming the Calzetti law for dust extinction, solar metallic- 
ity, and an underlying stellar population with 300 Myr of 
constant star-formation) to preclude possible systematic 
adjustment errors. The resulting Eb-v values confirm the 
trend indicated by Fig. ED 
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Fig. 18. The dependence of the strength of the contin- 
uum slope parameter j3 (see Leitherer et al. 2002J on the 
Lya equivalent width. Symbols are as in Fig. 1161 The (3 
uncertainties were obtained from the scatter of the values 
derived for the individual spectra. The values given in S03 
were calculated from the colour excess Eb-v, whereas the 
FDF values were derived by fitting the composite spectra. 
Hence, the S03 values (crosses) were corrected assuming 
a constant offset. The shift (A/3 = -0.92 ± 0.05) was es- 
timated from the (3 measured in the published composite 
spectrum of the total S03 sample. 

According to IS03I the (3 — WLya relation can be ex- 
plained assuming that starburst regions are covered by 
dusty gas clouds, which cause a decrease of the escape 
probability of Lya photons (see above) and an increase of 
the extinction of the stellar continua by dust. To allow sig- 
nificant Lya emission from galaxies at z < 3, the increase 
of the dust content in luminous galaxies at lower redshifts 
must not lead to a strong decrease of the escape probabil- 
ity of Lya photons. Hence, the observed evolution of the 
(3 — WLya relation requires that the amount of neutral gas 
grows more slowly than the dust content with decreasing 
redshift and/or the typical geometric distribution and the 
velocity field of the interstellar medium have to change 
suitably. 

The nebular emission feature C III] AA 1907, 1909 
traces the electron temperature in HII regions, which de- 
pends on the composition of the stellar population and 
the metallicity (Heckman et al. I1998|l . Fig. 1191 shows that 
the FDF galaxies (open circles) follow well the relation 
between Lya and C III] observed in the S03 sample, indi- 
cating strong C III] emission, and thus hot H II regions in 
the case of Lya-bright galaxies. 

Finally, we investigated the dependence of the equiv- 
alent width of the He II A 1640 emission feature on WL yQ - 



Fig. 19. The dependence of the strength of 
C III] AA 1907, 1909 on the Lya equivalent width. 
The open circles represent the values measured in the 
FDF composite spectra for different Lya strengths. For 
comparison the corresponding data from S03 are also 
included (crosses). 

Hell A 1640 can form only if a hard continuum is present. 
Excluding AGNs, significant Hell emission is expected 
in the presence of Wolf-Rayet stars only (Schaerer 2003 ) . 
Since Wolf-Rayet stars require a young starburst (age in- 
terval 3 < t < 7Myr for solar metallicity), Hell is an 
important indicator of the recent star-formation history 
and the metallicity. Fig.l20lshows that for the majority of 
FDF spectra the EW of Hell is close to -1 A. This re- 
sult agrees well with Wh c ii derived for the S03 composite 
spectrum for the whole sample (cross). On the other hand, 
an inspection of the individual FDF spectra show partic- 
ularly large variations of the Hell feature. As an exam- 
ple we note that FDF-5903 (Fig.[SJc)) has a significantly 
stronger Hell emission (Wkeii = — 1.6 ± 0.2 A) than the 
average while other galaxies show He II emission statisti- 
cally significantly below the average. Furthermore, Fig. 1201 
suggests that strong Lya emitters have particularly high 
WhoIi- However, because of the large error in WhoII of the 
Whya = —53 A point, this result is not statistically signif- 
icant. These observations may indicate that, like in the 
local universe, there are high-redshift 'Wolf-Rayet galax- 
ies' having exceptionally young starburst populations. 

9.6. Relations with the continuum slope (3 

In principle, the UV continuum slope [3 can be affected by 
various effects, such as the star-formation history (which 
determines the composition of the stellar population) , the 
initial mass function, the chemical composition, dust ex- 
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Fig. 20. The dependence of the strength of He II A 1640 on 
the Lya equivalent width. The open circles represent the 
values measured in the FDF composite spectra for differ- 
ent Lya strengths. The cross indicates the He II equivalent 
width of the composite spectrum of the S03 sample. 



tinction and the contribution of a possible AGN. However, 
for starburst galaxies showing a very young stellar popula- 
tion the predominant effect is expected to be dust redden- 
ing (see Calzetti et al. 11551 Leitherer et al. 115351 12UU2|) . 
This implies that [3 depends on the distribution of dust 
in the galaxy and on the metallicity (see Heckman et al. 
1998 ), which has an important impact on the dust forma- 
tion. 

In Table 13/3 was found to be strongly redshift depen- 
dent, evolving from (3 — —1.8 at z ~ 3.3 to (3 = —0.6 at 
z ~ 2.4. To study this effect in more detail we also calcu- 
lated the (3 values of the individual galaxies. The resulting 
Fig. 1211 confirms the change of (3 with redshift. 

To exclude that these results are due to a luminosity- 
related selection effect, we plotted in Fig. 122/3 as a func- 
tion of the galaxy luminosity at 1500 A (see Fig. Ill|) . 
As shown by the figure, the z > 3 objects show a dis- 
tinctly higher average luminosity ((logLi5oo[W/A]) = 
34.35 ± 0.04 in comparison to 34.14 ± 0.03 for the lower 
redshift range). This effect is mainly due to our selec- 
tion criteria, causing an increase of the luminosity limit 
with increasing redshift. Thus, we cannot exclude the exis- 
tence of reddened faint z > 3 galaxies. On the other hand, 
the z < 3 sample contains only a few 'blue' galaxies, al- 
though the brightness limits would allow their detection. 
Furthermore, a different sample selection would not in- 
crease the number of such galaxies significantly, since only 
a few luminous FDF objects have suitable magnitudes and 
colours. Finally, the / magnitude based sample selection 
does not introduce a significant /3-related selection effect. 
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Fig. 21. UV slope (3 between 1200 and 1800 A (Leitherer 
et al. I2002|l as a function of redshift for the FDF galaxy 
subsample used in the calculation of the composite spectra 
(see Sect. 19 .2|) . Typical errors of (3 are around 0.1. 



Assuming no evolution, the ratio of 'blue' to 'red' galax- 
ies should be almost the same for z ~ 2.3 and z ~ 3.3, 
since I(z ~ 3.3) - I{z ~ 2.3) w /(1800A)//(2400A) 
is similar for both galaxy types (compare Type IV and 
V in Fig. [3J. Consequently, the conspicuous difference 
between the mean (3 of the logi 1500 [W/A] > 34.2 ob- 
jects (((3) = —0.58 ± 0.13 for z < 3 in comparison to 
(/?) = —1.74 ±0.13 for z > 3) is obviously real and cannot 
be explained by the selection procedure applied. 

The very different continuum slopes for similar abso- 
lute brightnesses suggest that there is a strong evolution 
of the intrinsic UV luminosities of the galaxies. Hence, we 
estimated this quantity following Leitherer et al. (j2002). 
Taking 



Aisoo = 2.19 08 - fa) 



(10) 



(A1500 being the attenuation at 1500 A in mag), /3 = —2.5 
and 



log Ltjv = log L 



1500 



OA A 



1500 



3.2, 



(11) 



we obtained the values plotted in Fig. As expected, 
there is a significant redshift dependence of the extinction- 
corrected total UV luminosity, leading to very luminous 
galaxies (£uv > 10 12 Lq) at z < 3, while such objects 
are rare at higher rcdshifts. The galaxies at z < 2.5 
with Ltjv ~ 10 13 Lq are particularly interesting, since 
they are as bright as the most luminous star-forming 
galaxies known, implying star-formation rates of several 
100 Mq/yt and more. Fig.[2Slcan be explained in the sense 
that the masses of bright starbursts grow with cosmic age. 
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Fig. 22. UV slope f3 as a function of the luminosities at 
1500 A for the FDF galaxy subsample used in the calcu- 
lation of the composite spectra. Squares mark galaxies at 
z < 3. Triangles indicate galaxies at z > 3. Typical errors 
of P are around 0.1. 

Next, we analysed the dependence of j3 on other spec- 
tral properties. For this, we divided our sample of 64 
galaxies into four equally-occupied (3 intervals and con- 
structed composite spectra. Further /3-dependent mean 
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Fig. 23. Total extinction-corrected UV luminosity in solar 
units as a function of redshift for the FDF galaxy subsam- 
ple used in the calculation of the composite spectra. 



Fig. 24. The dependence of the strength of low ionisation 
lines (see Fig. I16|) on the continuum slope (3. The open 
circles indicate values measured in FDF composite spec- 
tra representative of the different quartiles (containing 16 
objects each) of the f3 distribution derived from the indi- 
vidual spectra. The subsamples of the ranges 2 < z < 3 
(filled squares) and 3 < z < 4 (filled triangles) were di- 
vided into two parts of equal number and averaged using 
P as classification parameter. Data points representing the 
six steepest continua (f3 < —1.3) at z < 3 and the five flat- 
test continua (/3 > — 1.3) at z > 3 are marked by a small 
filled square and triangle, respectively. For comparison the 
corresponding data from S03 are also shown (crosses). 

spectra were calculated by halving the subsamples of the 
ranges 2 < z < 3 and 3 < z < 4 according to the con- 
tinuum slope. The strong (3 evolution made it necessary 
to construct composite spectra including the six steepest 
spectra (f3 < —1.3) with z < 3 and the five flattest ones 
(j3 > —1.3) with z > 3, respectively. 

In Fig. [2] we show the average EW of the six strong 
low-ionisation lines (see Fig. I16|) as a function of p. The 
figure indicates a monotonic relation, characterised by a 
decrease of VFlis with bluer continua. A slight redshift de- 
pendence of the relation cannot be excluded, since Wlis 
seems to be smaller at constant (3 towards lower redshifts. 
Heckman et al. (1998) conjecture that the Wlis — P rela- 
tion could be caused by a relation between dust extinction 
and turbulent velocity of the interstellar medium. This 
could be due to the dustier galaxies being more massive 
and showing stronger and more violent star-formation. 
According to SOS] the Wlis — P relation may be due to 
the change of the covering fraction of the hot stars by in- 
terstellar clouds consisting of neutral gas and dust, which 
gives a more direct link between both quantities than the 
former explanation. 
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Fig. 25. The dependence of the CIV equivalent width on 
the continuum slope parameter /3. The data points were 
derived from composite spectra for different (3 and z in- 
tervals. Symbols are as in Fig. 1241 

Fig. [23 indicates a relatively tight relation between 
Wciv and (3, which does not show a redshift evolution. 
In part the strong decrease of Wciv with (3 for j3 < —1.5 
can be due to a C IV emission component. Fig.l25lsupports 
the assumption that the dust reddening depends partly on 
the metallicity, which is traced by C IV (see Heckman et 
al.EUBI Leitherer et al. l2UTTTl Mehlert et al. l2UU2jl . 

10. Implications and conclusions 

As described in the preceding sections, in addition to pro- 
viding valuable information on the spectral properties of 
distant galaxies, the FDF spectroscopic survey also pro- 
vided interesting new data on the evolution of the prop- 
erties of bright starburst galaxies 2 < z < 5. The main 
results on this redshift evolution can be summarised as 
follows: 

— Although the spectra of the starburst galaxies are ba- 
sically very similar at all observed redshifts, there is a 
tendency for the net line absorption to become weaker 
with increasing redshift. In particular, our galaxy spec- 
tra with z > 3 show on average weaker absorption 
and/or stronger emission components. The effect is 
particularly conspicuous for the Lya line and the C IV 
resonance doublet. 

— The UV continuum slope tends to become flatter with 
decreasing redshift. Since for bright starburst galaxies 
this slope is mainly determined by the internal redden- 
ing, this correlation also indicates an increase of dust 
reddening with cosmic age. 



Our data show an anticorrelation of low-ionisation 
interstellar absorption lines with the Lya emission 
strength (reported already bv lS03|l and a positive cor- 
relation of these lines with the reddening. These corre- 
lations do not show a significant redshift dependence. 
The C IV absorption strength (while clearly dependent 
on the redshift) shows only little correlation with the 
Lya emission or absorption strength. On the other 
hand, a positive, redshift-independent correlation be- 
tween the C IV absorption strength and the dust red- 
dening was found. 

The anticorrelation between the reddening and the 
Lya emission strength shows a striking redshift de- 
pendence, at least for absorption-dominated Lya pro- 
files, where the average UV continuum flattens with 
decreasing redshift. 

Although our flux-limited spectroscopic sample is ob- 
servationally biased towards absolutely brighter ob- 
jects at higher redshifts, the average intrinsic UV lumi- 
nosity of the starburst galaxies in our sample decreases 
with redshift for 2 < z < 4. 



As shown by Mehlert et al. (2002) the increase of the 
C IV absorption strength with decreasing redshift can be 
readily explained by the cosmic heavy element enrichment 
history of starburst galaxies in an evolving universe. The 
cosmic chemical evolution also provides a plausible expla- 
nation for the observed dependence of the UV continuum 
slope (3 and the Lya emission strength on the redshift. 
For a lower metallicity we expect a lower dust content (and 
consequently a lower reddening) , a higher escape probabil- 
ity for Lya photons, and higher temperatures of the radia- 
tion field and gas of the starburst galaxies. This prediction 
is in good agreement with the observed steeper continuum 
slope and (on average) larger Lya emission observed for 
z > 3. 

On the other hand, the Lya emission shows a rela- 
tively strong correlation with the low-ionisation absorp- 
tion lines as compared to the relatively poor correlation 
with the C IV absorption. Since Wlis traces the properties 
of (dusty) interstellar HI clouds fsee lS03|l . a strong influ- 
ence of the velocity field and/or the geometric structure 
of the interstellar medium (and resulting covering frac- 
tions) on the observed Lya flux and its evolution can be 
assumed. The observational result that a significant frac- 
tion of the galaxies with z > 3 show Lya in absorption, 
although the average reddening is distinctly lower than for 
the corresponding 2 < z < 3 galaxies, could also be ex- 
plained by an evolution of the properties of the interstellar 
medium. 

The decrease of the average Lya emission and the con- 
spicuous enhancement of the dust reddening with decreas- 
ing redshift may be partly due to an increase of the average 
cold gas and dust mass of the galaxies, which impedes the 
escape of UV photons from the galaxy by an enhanced ob- 
scuration of the luminous hot stars and H II regions in the 
starburst cores. Assuming that the intrinsic UV luminosi- 
ties of the high-redshift galaxies are characteristic for their 
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starburst mass, the observed increase of the mean abso- 
lute luminosity with decreasing redshift of the galaxies in 
our 2 < z < 4 sample suggests mass evolution as well. 
This result can be interpreted within current hierarchical 
models for the formation and evolution of galaxies. 

Summarising, wc conclude that the observed evolution 
with redshift of the basic properties of the galaxies in the 
FDF spectroscopic survey sample can be well explained 
on the basis of present theoretical concepts and theories of 
the formation and chemical evolution of galaxies. Hence, 
the data described in this paper provide strong further 
support for these theoretical concepts and models as well 
as a potential basis for additional investigations of the 
details of such models. 
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